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Microscopic investigation is performed for intersite multipolar interactions in the orbitally 
degenerate Anderson lattice, with CeBg taken as an exemplary target. In addition to the 
/° intermediate state, p Hund's-rule ground states are included as intermediate states for 
the interactions. The conduction-band states are taken as plane waves and the hybridization 
as spherically symmetric. The spatial dependences of multipolar interactions are given by 
the relative weight of partial wave components along the pair of sites. It is clarified how the 
the anisotropy arises in the interactions depending on the orbital degeneracy and the spatial 
configuration. The stability of the antiferro-quadrupole order in the phase II of CeBg 
is consistent with our model. Moreover, the pseudo-dipole interactions follow a tendency 
required by the phenomenological model for the phase III. 

KEYWORDS: quadrupole order, orbital degeneracy, pseudo-dipole, octupole, Anderson lattice, 
RKKY interaction 

1. Introduction 

Rare-earth compounds with orbital degrees of freedom have been attracting much interest. 
Among them, CeBg shows antiferro-quadrupolar (AFQ) ordering below Tq = 3.2K called the 
phase II, and antiferro (AF) magnetic ordering below T/v = 2.3K called the phase III.^'^^ 
This AFQ ordering has been identified as of quadrupoles of Fs type.^^ On the other hand, 
the phase III has a complicated magnetic structure with double fc, i.e., k\ = (1/4,1/4,1/2) 
and k\* = (1/4, —1/4, 1/2) in units of 27r/a with a being lattice constant, as has been probed 
by neutron scattering.^) The magnetic structure of the phase III is such that parallel and 
anti-parallel magnetic moments of next-nearest-neighbor pairs are equally populated. Thus 
the magnetic structure has no gain of energy from the isotropic exchange interaction acting 
on nearest neighbors. In order to explain the stability of the structure, pseudo-dipole-type 
couplings have been proposed between the next-nearest-neighbors. However, microscopic 
justification of this type of anisotropic interactions is lacking. 

There are two different sources of the intersite exchange interaction: the on-site Coulomb 
exchange interaction mainly between 4/ and 5d electrons in rare-earth systems, and hybridiza- 
tion between 4/ and ligand electrons. The Coulomb exchange always favors on-site parallel 
spins. Early studies derived its consequences neglecting the orbital degeneracy, and assuming 
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the spherically symmetric conduction states.^) The intersite quadrupolar interaction was de- 
rived later for the model with the Coulomb interaction as the source. The most elaborate 
contribution for the Coulomb model is made by ref.7, which includes the multipoles of 4/ elec- 
trons and the result of band structure calculation. On the other hand, intersite interactions 
originating from hybridization have been considered in much less detail. Theoretical litera- 
ture often concerns with isotropic interactions only, which result from neglecting the multiplet 
splitting in the configurations.^'^) The intersite interaction from band structure calculation 
was derived in ref.ll. Recently Shiba et al. derived the effective interionic Hamiltonian with 
the full use of cubic symmetry of the conduction states, but only the intermediate state 
was included. 

It has been shown that the anisotropic spin exchange arises by cooperation of the spin- 
orbit coupling and the Coulomb exchange interaction.^' However, the anisotropic exchange 
from hybridization has hardly been studied. Motivated by this situation, we investigate in 
this paper the multipolar interactions induced by hybridization microscopically. We include 
the spin-orbit interaction in the intermediate states to derive possible anisotropy. Specifically 
we introduce the Anderson-type lattice model with orbital degeneracy, in which not only 
but also Hund's-rule ground states are considered as intermediate states. 

In order to clarify the origin of the anisotropy and the difference between multipolar inter- 
actions in a simple manner, we take a model which has spherically symmetric hybridization 
around each Ce ion, and which takes the wave functions of the conduction band as plane 
waves. For the ground state of each Af^ configuration, we take the Fs crystalline electric field 
(CEF) state, bearing CeBp in mind. Then we obtain the effective interionic Hamiltonian by 
4th order perturbation with respect to hybridization. The high symmetry around each Ce ion 
allows us to simplify the analysis significantly. 

As is well known, the intermediate state tends to align the spins of / and conduc- 
tion electrons anti-parallel. On the other hand, the Hund's-rule correlation works in the 
intermediate states. The lowest intermediate multiplet, referred to as hereafter, tends to 
align the and conduction spins parallel because the total spin 5 = 1 is realized in the 
configuration. This competition causes decrease of the spin exchange as will be investigated 
in detail later. 

This paper is organized as follows. In §2 we introduce the model modifying the periodic 
Anderson Hamiltonian with orbital degeneracy, and calculate the effective interionic Hamil- 
tonian by the 4th order perturbation. In §3 we evaluate the multipolar coupling constants 
numerically, and analyze them in detail. We comment on possible relevance of the results to 
actual systems. The summary is given in the final section. Appendix derives analytically the 
range functions of the intersite interaction within our model. 
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Each Ce^+ ion has one 4/ electron, and the spin-orbit coupling leads to the total angular 
momentum ji = 5/2 in the ground state of the ion. The operator fma{i) creates an 4/ electron 
at site Hi with the z-component m of the orbital angular momentum Z = 3 and spin cr. We 
take conduction-band wave functions as plane waves, and introduce the creation operator 
with wave number k and spin a. Then our model is given by 



where 



H = Hc + Hf + ifhyb + -f^Coulomb, (1) 



Hc=^ekcl^Cka, (2) 



k.a 



A 



Hf = Yl EffLaii)fma{i) + - J] Sm\L\r, 3m) ■ f}n„{i)f,n'a'{i), (3) 

i,m,<7 i,m,m' ,cr,(7' 



hyb=V^ E [^fe>^3*mWe^'=-'^/L«Cfea + H.c.]. (4) 

i,fc,m,(T 

Here A is the spin-orbit coupling constant, while L and a denote the orbital angular momen- 
tum operator and Pauli matrices respectively. The solid angle of k is represented by Q.^. We 
assume that combination of Ef and -ffcouiomb is such that the 4/-'^ state is stable against 4/° 
and 4/^ states. 

As the intermediate states with configuration, we consider only the Hund's-rule ground 
state /h in order to emphasize the origin of the anisotropy. In general, is characterized by 
the quantum numbers L, S, J and Jz = M. We take L = 5, S = 1 and J = 4 according to the 
Hund's-rule. Then we replace Hf + i?Couiomb as follows: 

Hf + i^Coulomb 

^ ^ £;o |z; 0) (z; 0| + ^ £^1 |i; iimi) {i; jimi | + ^ ^2 K; {LS)JM) {i; {LS)JM\ , (5) 

i i,mi i,M 

where Ei is the energy of the 4/^ ground state, while Eq and E2 are the energies of and 
f^, respectively. Experimentally, the position of the 4/^ level is about 3 to 4 eV below the 
Fermi surface in CeBg as probed by photoemission.^'^'^"^) On the other hand, the 4/^ levels 
have been probed by Bremsstrahlung, and they are also 3 to 4 eV above the Fermi level. 
Thus naively thinking, intermediate states 4/*^ and 4/^ have comparable contribution to the 
effective exchange interaction. It is not clear, however, the weight of the 4/^ relative to other 
contributions. Our theory with only 4/^ taken into account may give the upper bound of the 
anisotropy in the exchange. 

In eq.(4), we then need to consider the matrix elements between /°, and 4/^ states: 

-f^^hyb -ffhyb(/°) + -f^hyb(/^), (6) 
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where i?hyb(/'') is the hybridization Hamihonian between Af^ state and f^, while -ffhybl/n) 
is that between 4/^ state and /y. They are represented by 

fin 



i,k mi,m,(7,ai 



X {jimi\3mi-ai,^ai){mi-ai,ai\fl^j0) 

X |i; jiTOi)(i;0|cfe^ + H.C.] , (7) 
^hyb(/^) = y^E E [^fc^^mWe*^ 

i,k M,m\,m,a,a\ 

X ((L,S')JM|/^^|mi-(Ti,ai)(3mi-(Ji, ^ ai\jim,i) 

X \i;{L,S)JM){i-jimi\cua + ^-C-]- (8) 

^.i Effective on-site exchange Hamiltonian 

Using eqs.(7) and (8), the effective on-site exchange Hamiltonian is obtained by the second 
order perturbation theory as:^^) 

Hexch = Hexchif^) + -f^exchC/n)) 

gexch(/°) = -| E E -f'-e ''^"'"''"^ E nmA^>^)ySms{^k') 
i,k,k' ° ^ All m's,a's 

X (mi-c7i,cri|/44(,JO)(0|/^3^3|m2-(T2,a2) 
X 0'imi|3mi-(7i, ^ (7i)(3m2-(T2, ^ cr2|iim2) 

X \i;jimi){i;jim2\cl,^_^Cka4, (9) 

i,fe,fe' ^ ^ All m's,(T's 



^(mi-ai,cTi|/^3,3|(LS)JM)((L5)JM|/t,^,Jm2-a2,c72) 



X 
M 



X 0imi|3mi-(7i, ^ (7i)(3m2-(T2, ^ cr2|iim2) 

X jim2|c^,^gCfeo-4, (10) 

with being the Fermi energy. In the above, -ffexch(/'^) and i?exch(/H) the effective 
exchange Hamiltonian with the virtual states /° and /g, respectively. 

Let us inspect the sign of the effective Hamiltonian (9) and (10). If the multiplet splitting 
is neglected in the 4/^ states, both /° and 4/^ give the antiferro-magnetic spin exchange 
interaction, as is well known for the Anderson model. If the Hund's-rule splitting is included, 
/h favors the ferromagnetic spin exchange. This is demonstrated as follows: from the effective 
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Table I. Signs of exchange interactions between conduction and 4/ electrons with given intermediate 
states where F and AF mean ferro- and antiferro-couplings, respectively. The spin triplet is in the 
column f^{S = 1), and the singlet in = 0). The multiplct is a special case of f'^{S = 1). 



exchange 


/° 


f\S = 1) 


f{S = 0) 


spin 


AF 


F 


AF 


orbit 


AF 


AF 


F 


spin and orbit 


AF 


AF 


AF 



Hamiltonian (10), spin exchange corresponds to the case of ui = o"4, cr2 = o'3, mi—ai = m2— o"2(= 
m^) and rris = m^. In this case, the relevant matrix elements in eq.(lO) combine to give 

^{m,ai\fm,as\iLS)JM){iLS)JM\fl^,^\m5a3) 



M 



1[1 _ (_l)3+3-L+i + i-5]2(_^^l + i-5 



M 



{JM\LML,SMs){LML\3ms,3m5){SMs\^as, ^a^] 



(11) 



Ml, Ms 

where we have used the following property of the Clebsch-Gordan coefficients: 



(j2m2,jimi|j3m3) = (-iyi+^2 ^3(jimi, j2m2|j3?T^3)- (12) 

It is seen that L + S must be even for eq.(ll) to be nonzero. Since the Clebsch-Gordan 
coefficients are real, the right hand side of eq.(ll) becomes positive for 5 = 1, which is a 
property of f^. Considering that the energy denominator in eq.(lO) is negative, we derive a 
ferromagnetic on-site exchange, which is common with the double exchange mechanism. 

In the case of the orbital exchange, we have a result similar to eq.(ll) except that the sign 
factor is replaced by (— 1)^+^"^. With 5" = 1, wc obtain the negative sign factor, which brings 
about an antiferro-type interaction for the orbital exchange. The case for the simultaneous 
exchange of spin and orbit has the sign factor (— l)3+3--f'+2 + 2~'^ which is also negative for 
Thus the corresponding interaction is again of antiferro-type. The 4/^ multiplets have excited 
states with S = 0. With S" = in the intermediate states we obtain different signs depending 
on types of the exchange. Table I summarizes the signs for various intermediate states. We 
note that these consequences come from the antisymmetry of the 4/^ wave function. 

To show the symmetry in eqs.(9) and (10) more clearly, we introduce the creation operator 
of the total angular momentum j as 

fL = J2fl<r{^^'^r\j^)- (13) 

m',cr 

Using the Wigner-Eckart theorem, the matrix element of /j^ is rewritten by its reduced matrix 
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element as 

ihmilfjjO) = (-l)^--i ( ^ l) Oi||/j||0). (14) 

y— mi m UJ 

Using this relation and introducing another angular momentum Pi and the 6-j symbol/^"^^^ 
we write the part including the matrix elements of fermion operators as 

O'l ^1 1 4m J 0) (0 1 fjumu \h 

= -{-l)n-rm+h-m, y2 A{f;jt,ju,Pi)[Pi]( A JiUjt Pi (15) 

p^ Q^ \-mi Qi 7712/ V~"^* Qi "^"Z 



J2{hrni\fj„rnJiLS)JM){iLS)JM\fl^Jnrn2) 



M 

p^Q^ \-mi Ql 7772/ V""^* Qi "^V 

where [j] = 2j + 1 and 

A(/°;jt,i„,Pi)=-Oi||4l|0)(ii||4||0)*(-ir+o+^-^J ^M, (17) 

[jt Ju J 

A{fl-JtJu,Pi)={{LS)J\\fl\\h)^iLS)J\\f^^^^^^^^ ^l]. (18) 

[jt Ju J ) 

The curly brackets in eqs.(17) and (18) are the 6-j symbols. We note that A(^; j, j' , P) includes 
information of wave functions of the ground state and the excited state specified by ^. The 
transformation described above is best handled by a diagrammatic technique for addition of 
angular momenta. Using the theory of angular momentum and the antisymmetry of wave 
functions, we find that 

(/ll/jl|0) = -<5,',,b-r/', (19) 

'3 3 L 

((L5)J||4||j') = (-l)^'-^"+-'^/2[L,5,i,/,J]V2<ji/2 1/2 s\, (20) 

j f 

where the curly bracket in cq.(20) is the 9-j symbol, ^^"^^^ and the shorthand notation 
bi> ■ ■ ■ > in] = (2ji + 1) • • • {2jn + 1) is used. In this way we obtain the following form: 

5/2,7/2 2n f"Jl . . p. 

^^„. = >,i=i: Y. Y. Y. i^ /'^f ' I') 

i,fc,fc' 3UU Pi=0 excvsi; 



X 

a.a' mt ,mu 



X (-l)^'^-™^+^'*-"^*(jt777t|3777t - a, ^a){3mu - a', \cj'\jumu) 
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X \i;jimi){i;jim2\c]^,^,Cka, (21) 

where is the index of intermediate states at Ri, which is either /° or f^. The excitation 
energy E'exc is defined as follows: 

E,^,{f) = Ei-Eo-eF, (22) 

E^^,{fl) = Ei-E2 + eF. (23) 

In place of the set of (2ji + 1)^ operators \jimi){jim2\ with mi,m2 = —ji, ■ ■ ■ ,ji, we can 
use more conveniently the set of irreducible tensor operators Tq, with Q = —P, ... ,P, and 
P = 0, . . . , 2ji.^'^^) The irreducible tensor operator is related to \i; jimi) {i; jim2\ by,^) 

\i;nrm){i;nm2\ = {-1^-- ^^j, ( ^^(0, (24) 

P^QQ^_P {jiW^ \\n) \-mi Q 1712/ 

where (P = 1, . . . ,2ji) is a reduced matrix element. 

By inserting eq.(24) into eq.(21), we obtain another expression for the on-site exchange 
Hamiltonian 

5/2,7/2 2n fJl ..^ . . p. 

^exch = 2^ 2^ 



" ' i?. IT ^-te'"'"^''"*) 

ct,<t' mt,mu 

X {jtmt\2,mt - a, ■^C7)(3m„ - a', -^o-'ljumu) 



2 " ' 2 

Pi 

X 

where the orthogonality relation of 3-j symbols is used.^^^ 



E (-ir-M '^VSW^V^'^^^' (25) 



S.S Effective intersite interaction 



In order to derive the intersite interaction we work with the waver-number representation 
of the on-site interaction: 

i,ii,Pi,Qi k,k',a,(j' 

where the matrix element ^k^^k^i^i) is determined by comparison with eq.(25). Then the 
effective intersite interaction between two Ce ions, which are located at Ri and R2, is given 
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by the second-order perturbation theory with respect to H^xch as 

k,k',a,a' ^ ^ Pl,P2,Ql,Q2 

xE<::?^.(^i)<^?v(6)- (27) 

a, a 

We replace the sum over k by integrals, 

E ^ ^ (1^) J ''^dk J dOfe. (28) 
The angular integral is of the form 



!■ 



d^kYiruf {^k)Y^mi {^k) exp(±ifc • R), (29) 
for which we use the partial- wave expansion of a plane wave: 

exp(±zfe • r) = 'iT:Y,{^i)'jl{kr)Y,YL{^k)Yim{r), (30) 

I m 

where f = r/r denotes the unit vector in the direction of r. Furthermore the decomposition 
of the product of spherical harmonics with the same argument:^^) 



0> 



X 

M 

Then we obtain 



1 °° /l 
(i^)ll,^,(i^) = ^ j;(-l)^[L,Zi,Z2]^ ( 

E t '0>^.M(n). (31) 



/ 



dnkYimf{^k)Y3m,{^k) expi±ik ■ R) 



1=0 VO 0/ m=-l V™/ 

In eq.(32) the first 3-j symbol vanishes for odd l,^^^ hence we have to consider even I only. 

In the next stage we take summation over intermediate angular momenta, and simplify 
the results by using symmetry properties. For example, the summation over a, a' in eq.(25) 
gives rise to a 6-j symbol:^^'^^) 

h 3 \ 



yZijtmtlSmt - a, Ja)(3m, - a, la|j,m,)(-l)3+(— ) ( ^ 



— a m —{ms — (j)l 

{-rr'^'^hjs,k]'''\^ ' \ {-ly^-^ ( V (33) 

t jt js 2 ) \mt m -rus J 
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and the summation over internal azimuthal quantum numbers yields the 9-j symbol:^^'"^^^ 

jt Pi ju \ 



^_l'j{jt-mt)+{js-ms)+{ju-mu)+{jv-mv)+{L-M) 
rus ,mt ,mu ,mv ,m,m' 

jt h js 



{-mt Qi mj 



X 



jv I: 



2 Ju 

I 



\mt m —TxisJ \rny m — m„y 



r 3v P2 js] 

[-my Q2 ms) 



—M m m'l 



(-1) 



P2+L+h+jt-3v 



(Pi P2 L\ 
[Qi Q2 Ml 



(34) 



Pi ju jt\ 
' P2 jv js 

L h 

We then take the integral over radial momenta k,k' . This is given by 

QkhikpR) = (kpRr'^ x^dx I ^ ^ [ji^ {x)ji^ {x') + ji-, {x')ji^ {x) } , (35) 

where ji{x) is the spherical Bcssel function of order I and the energy dispersion for the free 
electron = h?k'^/2m is inserted. Note that Qi^i2{kFR) is invariant under interchange of li 
and I2 by its definition. Appendix gives the analytic results for Qi^i2{kFR) for all relevant 
values of lil2- Putting all together we finally obtain the following form: 

2ii 



5 1 

2 ' 2 



= E E E 

47r2 ^^^^2 h,l2,L=0 Pl,P2=0js,jt,ju,jv 

^ H^i;jt,ju,Pi)H^2;jv,js,P2) (I, T> ^ 

-C'exc(,?lj-C'exc(,?2j 
x(-l)('l+'2)/2-i.+.*-i„-i.-Pl+l[3]2[^^^;^^p^^p^][_^.^^_^.^^_^.^^_^-^^^]i 



(3 h 3\ (3 h 3\ L h M 1 3 3 Zil [3 3 h 



VP 0/ \p 0/ \0 Oy 

{TP^{1)TP^{2)Yl{R^2)}1 



\jt Js 2 ) -^^ 2 > 



Pi 


Ju 


jt 


P2 


jv 


js 


L 


h 


h 



(36) 



(ii||T^MIji)(jil|r^^l|ji) ' 
where R12 = Ri — R2- In eq.(36) we have used the fact that L,li,l2 are even integers, and 
introduced Ea = h^/2ma'^. In eq.(36), {• • • }o denotes a tensor component of rank zero, i.e., a 
scalar,^' It can also be written as 

{T^^{1)T^^2)Yl{Ri2)}^o 



(_i)-Pi+-P2+i 

QuQ2,M 



T^imSimYLM{Ri2). 



(37) 



(Pi P2 L\ 

VQi Q2 m) 

Among various components of the Hamiltonian, those terms with L = represents the 
isotropic interactions. Other terms with L ^ represents anisotropic interactions which de- 
pend on the relative angular position of the two ions. For example, the pseudo-dipole in- 
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teraction corresponds to the case Pi = P2 = 1 and L = 2.^^ Let us choose i?i2 along the 
quantization axis z. Then we have R12 = z, and Ylm (-R12) vanishes except for M = 0. Then 
the dipole interaction has the fohowing form 

K{Ru)[3j,il)U2)-j{l)-ji2)]. (38) 

With K{Ri2) < 0, this interaction favors parallel moments along the z axis. For moments 
directed perpendicular to the z axis, the first term of eq.(38) vanishes and the second term 
favors antiparallel moments. This preference is just as in the case of real dipoles. On the other 
hand, with K{Ri2) > 0, moments along z tends to be antiparallel, and those perpendicular 
to z tends to be parallel. We refer to the case of K{Ri2) > as the "anti-dipole" interaction. 



2.3 Projection to crystal-field ground states 



So far we have not included any element of the crystal structure. In CeBg, the cubic 
symmetry around each trivalent Ce ion splits the sixfold degenerate ground state with ji = 5/2 
into a quartet called Fg, and a doublet called The crystalline electric field (CEF) splitting 
between Fg and F7 is about 540 K.^^) Hence at temperatures (< lOK) of our interest we can 
safely neglect the population in the F7 excited states. In considering the intersite interactions, 
we therefore pick out such component of 4/ wave functions that belong to the Fg states. In 
other words we perform projection of \jim) to Fg wave functions [va) (v = +,— and a =], |) 

by 



\jim) |w)(i/cr|m). 



(39) 



where specifies the orbital index and a the Kramers index, i.e., one of the time-reversal 
partners. We have omitted writing ji = 5/2 in the kets \jim) of the right-hand side. Explicitly 
the Fg quartet has the following wave functions: 



l+T) 
l+i) 
l-T) 
l-i) 



5 

+ 2 



+ 



+ 



+ - 



(40) 
(41) 
(42) 
(43) 



where each coefficient in front of |m) corresponds to {m\i'a){= (valm)). 

The transitions within the Fg have 15(= 4x4 — 1) components, and the operators de- 
scribing these transitions correspond to 15 multipole operators. Following the literature,^' 
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Table II. The list of multipole operators in the Fg subspace. 



r 


symmetry 


XT 


2u 


y/lExyz 




35 


(3^2 - r^)/2 






V3(x2 - y2)/2 




Aul 


X 






y 






z 




Au2 


x(5x2 - 3r2)/2 






y(5y2 _ 3r2)/2 






z(5z2--3r2)/2 




5u 


^/T5a;(y2 - •z^)/2 






Vl5yiz'^ - x2)/2 








c 




-\/3y2; 












-\/3x?/ 





we introduce two pseudo-spins cr and r in terms of the vector p^^/ of the Pauli matrices: 

^= kf^)Pr,r'(^'f^l> (44) 

a,T,T' 

^= (45) 

<t,<t',t 

Using these pseudo-spins, we can express a physical operator XT adapted to the point group 
symmetry,^' Here 7 specifies a component in the irreducible representation T with possible 
multiplicity. For completeness the multipole operators are summarized in Table II, where 
/ji,r],C are defined by 

fj, = (/x^, M^ M^) = (tV^, r^c7^ r^'cT^), (46) 
rj = [rf^rj^r]') = ^(-rV^+ VSr^a^, -rV^- VSrV, 2rV^), (47) 

^ = (C^C^C^) = ^(-V3TV^-rV^, V3TV2'-TV^2r^a^). (48) 

The subscript u represents the odd property under the time reversal, and g the even one. 
The operators belonging to V^g and V^g are classified into quadrupole operators. Under the 
cubic symmetry the dipole operator J are mixed with octupole operators. Hence cr and 77 have 
both dipole and octupole characters. The remaining representations T2u and V^u describe pure 
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Table III. Projection from into . 



'105 

1 B\/TD5(V+7a-) 



2 


±2 






2 


±1 






2 









3 


±3 


^ [3V3(±r + ^C^) + 7(Tr/" + ^r/^: 


) + (^C7^ + iay)] 


3 


±2 






3 


±1 


^ [15(±C - + 7v/3(±r/- + irf) 


+ V3(±c7^ + ia?')] 


3 





(7^^ + 





4 


±4 




4 


±3 




4 


±2 


i^(±2zM^ + r-) 


4 


±1 




4 








5 ±5 _^[V3(±^--i(i/) + (±^- + i^i/) + 2(Ta--zc7^)] 

5 ±4 + 

5 ±3 ^/J [3 V3(±C + K^) + 13(Tr/^ + ir?'^) + 2(±(7^ - )] 



5 =^2 lyh^"" 

5 ±1 ^^/5^[21(^C'^ + ^e) + ^/3(±??^ + ^r/^') + 14^/3(^a^-^^7^^)] 



36V77 



octupole operators. In Table III, the relations between irreducible tensor operators Tq and 
physical operators X'^ are summarized. As seen in the last section, the effective Hamiltonian 
(36) involves the product of operators with different ranks for each site. Hence after projection 
of the Hamiltonian into the Fg states, there should be couplings between different irreducible 
representations. Symmetry analysis of such couplings has been performed recently by Sakai 
et al. in analogy to the Slater-Koster scheme for the energy band theory.^^) In our case, the 
general discussion involves very many terms and is not illuminating. Therefore we focus in 
this paper only on the couplings within each irreducible representation T. 

In analogy to the dipole interaction given by eq.(38), we refer to the interactions of the 
form: 

^pseudo = K{Ri2){^[v{Ri) ■ Ri2][v{R2) ■ R12] - v{Ri) ■ V{R2)}, (49) 
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as the pseudo-dipolc interaction in the Tg subspacc. Here v is either of a, rj, C in Table II. 
The case K{Ri2) > is caUed the "anti-pseudo-dipole" interaction. 

3. Explicit Results for Multipolar Interactions 

In this section we derive strength and spatial dependence of multipolar interaction for a 
pair of Ce ions separated by R. We assume that the Fermi surface consists of a sphere with 
radius kp- Assuming the same number of electrons as that of the lattice sites with the spacing 
a, we obtain kpa = (Svr^)^/'^ ~ 3.09. Let represent the distance between the nearest- 
neighbor pair, and Rnnn between the next-nearest-neighbor pair. Then we have kpRnn ~ 3.09 
and kpRnnn ~ 4.37 for the simple cubic lattice. We note that a = A.IAA in CeBg. For 
evaluation of quadrupolar couplings, we take R along the (001) direction because couplings 
between nearest neighbors are important for AFQ ordering in CeBg. When R is parallel to 
the z-axis, Ylm{R) vanishes unless M = in eq.(37). This gives a high symmetry among 
components of multipole interactions. On the other hand, Ylm{R) with M ^ does not 
vanish for the directions of next-nearest neighbors. This gives a chance of larger anisotropy 
to the next-ncarcst-neighbor interaction. For the rest of this section, we define three cases 
depending on intermediate states on two sites: (a) Af^ intermediate states on both sites; (b) 
4/^ intermediate states on one site and 4/*^ on the other site; (c) 4/^ intermediate states on 
both sites. 

We comment on the relationship between the present model and the actual electronic 
structure of CcBg. The real Fermi surface of CeBg consists of three equivalent sheets, which are 
nearly spherical. It has been shown that not only the polarization within each pocket centered 
at X points in the Brillouin zone, but also the interpocket polarization plays an important role 
in determining the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction.^^) Thus, taking a 
single Fermi surface our model incorporates both intra- and inter-pocket polarizations in a very 
rough manner. Thus in a sense a model with the single Fermi surface is better than considering 
only the intrapocket polarization of three Fermi surfaces. In relation to the magnetic structures 
observed in phases II and III, we are particularly interested in the coupling constant for nearest 
and next-nearest neighbors. 

3.1 Quadrupole-quadrupole interaction 

We begin with the quadrupolar interactions, paying attention to the strength and sign 
at the distance between the nearest-neighbor pair. Namely, we are interested in whether the 
simple hybridization model can explain the antiferro-quadrupole order of the Fs type. We 
define the quadrupole-quadrupole interaction Hamiltonian as follows: 
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X Ei^3t(6,6;i2K(i)r^(2) + ^ <(a,6;«))""(i)/(2)], 



(50) 




where and /i" have been defined in Table II, and and ^2 specify either 4/" or 4/^ as 
intermediate states at Ri and R2, respectively. The dimensionless quantities L'p^(^i, ^2; -R) 
with r being 3g or 5g are determined by eq.(36) together with the projection to Fg. We use 
the matrix notation to represent the set of components. The form of the intersite inter- 
action depends on the excitation energies -E'cxc(O) since the i?-dependence of Dy{$.i,$.2', R) 
varies depending on ^1 and ^2- Namely, if the Ce ion tends to become tetravalent, we obtain 
|£'cxc(/'^)| < \Ecxc{f'^)\ and the intermediate state dominates. On the other hand, in the 
case of comparable Ecxc(f^) and i?cxc(/^)) all -Dr(Ci5 ^2! R) also have comparable multiplying 
factor. The situation in CeBg seems to correspond to the latter. 

Fig. 1 shows the spatial dependence of the quadrupole couplings D^g and D^g for R \\ 
(001). With the intermediate state f^, there are only two independent coupling constants for 
any distance along (001). The four-fold symmetry around the pair is high enough to lead to this 
property in consistency with the argument of ref.lO. In the case of as intermediate states, 
there are four different coupling constants. The anisotropy among them reflects the difference 
between the z component, which is parallel to R, and x or y component perpendicular to 
R. For the nearest neighbors we obtain D'^^ {^1,^2] zRnn) > for all a = x,y,z and all 
combinations of ^1 and ^2- This positive sign favors the simple AFQ ordering as observed in 
CeBg, and is also consistent with the temperature dependence of the elastic constants. 

On the other hand, the spatial dependence of D^g does not depend much on the intermedi- 
ate state or At the nearest-neighbor distance kpRnn ~ 3.09, we obtain the anisotropic 
behavior: D^^ < and Dfg > 0. This result is interpreted by the property of the transfer 
integral; if the 4/ wave functions at the site Ri and Rj are stretched toward each other, the 
transfer integral from Ri to Rj becomes large and the energy is lowered. In other words, 
we should have a ferro-quadrupole coupling between 02{i) and 02(i) short distance. The 
relation D^g < D^g will be discussed from a more mathematical point of view in §3.4. Since 
the absolute magnitude of D-^g is similar to that of D^g, it is not possible to explain the 
stability of the phase II in CeBg with the D^g order only from the above result. However, the 
observed antiferro-type order is consistent with our simple theory for the -Dgp interaction. 

3.2 Dipole-dipole interaction 

For convenience, we regard the vector representations F = 4ul, 4u2 as representing dipoles, 
although they actually mix with octupoles under the cubic symmetry. The pure octupoles cor- 
responds to F = 2u, 5u and are studied later. The intersite interaction between two "dipoles" 
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separated by R is then given by 

= y 1 

Ea f-^ -E'exc(Cl)-Sexc(C2) 

x,y,z 

X i:[<i(6,6;i«K( V(2) + i?4t2(ei,6;i2)<(i)/(2)], (5i) 

In addition to the nearest-neighbor interaction, we are much interested in the next-nearest- 
neighbor interaction in this case. Especially, we shall examine whether the sign of the pseudo- 
dipolar interaction is consistent with the one required by the phenomenological model.^^ 

Fig.2 shows the spatial dependences of D^ui- The direction of R in the region kpR < 4 
is parallel to (001), and the position of a nearest-neighbor at kpR 

nn 3.09 is included. For 

kpR > 4 we take R parallel to (101), which includes a next-nearest-neighbor position with 
nn 4.37. In Fig. 2(a) we obtain isotropic D4ui(^f^ , f^', zR^, which is consistent with 
the theory of Shiba et al. who started from a different model. We note that the isotropy in 
our case is generated only after the projection to the Fg CEF states. Before projection, a 
finite component with L = 2 in eq.(37) gives anisotropic dipole interaction. On the other 
hand, the next-nearest-neighbor interaction D iu\{f^ ■, ] Rnnn) becomes anisotropic, but is 
still isotropic in the xz-plane. We shall discuss the origin of the anisotropy in detail later, and 
only remark here that the /° intermediate state leads to ferromagnetic interaction at both 

Rnn ^-nd Rnnri' 

With included as intermediate states, the r4ui interaction becomes anisotropic even for 
R along (001) as seen in (b) and (c). In contrast with the case of Fs^ where the z-component fa- 
vors the ferro-quadrupole configuration, the z-component -Dfyi tends to be antiferro-magnetic 
with as seen in (b) and (c). 

We have checked that D 4ui{ f ^, f'^; zR) tends to cancel Diuiif^ ^ f'^'izR) and 
-D4ui(/h) /h! short distance (not shown in Figure). This is because the /° intermediate 

state tends to align the spins of / and conduction electrons antiparallel, while intermediate 
states tend to align them parallel. At short distance the propagation of conduction electrons 
just connects this tendency on two sites. 

Fig.3 shows the spatial dependence of D4u2 for R along (001) and (101). There appear off- 
diagonal elements D%^2 = -^Iu2 ^ along (101) direction. By symmetry we obtain D^^2 = 
^4u2 ^4u2 — ^4u2- Then it is convenient to take the principal axes as ey = {x + z)/\/2, 
e± = {x — z)/\/2 and y. Along the pair axis the relevant eigenvalue is given by -D4^2 + -^4u2) 
and along e± by Df^2 ~ ^tu2- The interaction is anisotropic even along (001). This can be 
understood if one recalls that F4u2 involves the orbital flip, in contrast with F4ui. Hence the 
anisotropy occurs more likely than in F4„i along (001). At the nearest-neighbor distance, the 
anisotropy is opposite to that of r4„i . We expect that the sum reduces the anisotropy of each 
contribution at Rnn- For the next-nearest-neighbor interaction, the anisotropy in Fig.3(a) with 
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intermediate state is opposite from that in (b) and (c). In (b) and (c), Df^2 ~ ^4m2 shows 
antiferro-magnetic coupling and is larger than the other components in absolute value. This 
situation corresponds to K{Ri2) < in eq.(49) and is consistent with the phenomenological 
theory for the phase III of CeBg.^^ However, the result in (a) cannot be interpreted in terms of 
eq.(49), since they differ significantly in all three principal axes. As for the behaviour at short 
range, which is not seen in Fig. 3, the intermediate states and give additive contribution 
to Dau2- This is because the on-site orbital exchange and the simultaneous exchange of orbit 
and spin, which are partly represented by D\u2-i have the same sign for and intermediate 
states. 



3.3 Octupole-octupole interaction 



We now turn to the octupolar interaction. This interaction plays an important role in 
understanding the NMR and neutron scattering results in CeBg consistently,^^ and also the 
phase diagram in magnetic fields.^' ^^-^ The interaction is given by 

Ho = y 1 

x,y,z 

X [D2n(ei,6;i?)r^(i)Tn2) + Yl (6,6;i^)C(i)C''(2)]. (52) 

a/3 

Fig.4 shows the spatial dependences of the octupole-octupole interactions and D^^ 
for R along (001) and (101). At the nearest-neighbor distance, we obtain > irre- 
spective of the intermediate states, which means the anti-ferro octupole coupling. On the 
other hand, D5u(/°, /°; i2) has a different spatial dependence from that of D5u{f^, f^; R) 
or D5u{f^, f^; R). At the next-nearest-neighbor distance, D^^—D^^ is negative and is larger 
than the other components in absolute value for cases (a), (b) and (c). This means that Fs^ 
moments on next-nearest-neighbor pair tend to be parallel to each other and along e±, namely 
perpendicular to R. This feature of " anti-pseudo-dipole" interaction is also consistent with 
the phenomenological theory for the transition from III to IIP under magnetic field. 

The above results for various kinds of multipolar interactions show that the dipolar, 
quadrupolar and octupolar interactions have the same order of magnitude. Especially the 
octupolar interaction between F2u type moments is as strong as the quadrupolar interaction 
between T^g type moments. If we consider only /°, the results are consistent with the previous 
symmetry analysis. In addition, the Fs^-type octupolar interactions show the range depen- 
dence similar to that of the Fs^-type quadrupolar interactions for a given set of intermediate 
states. 
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In order to understand better the nature of the intersite interaction, we apply the the 
idea of "bond density" of Shiba et al. to the present model. We can trace difference of 
spatial dependences, and identify the independent components describing the interaction. For 
simplicity we consider only the 4/° intermediate state. The bond density is defined by 

b{M, M',k,R) 

= Y,{hM\3M-a,^a){3M'-a, ^a\jiM') J dnkYij^_^{nk)YsM'-a(.^k)e"' '^ 

00 

= ^biiM,M',R)ji{kR). (53) 
/=o 

In the above we have used eq. (32) to decompose the bond density into partial- wave components 
bi- Each component bi is projected to the Fg subspace as 

Bi{va,v'a\R)= ^ {va\M)bi{M,M' ,R){M'\v' a') . (54) 

M,M' 

For deriving the spatial dependence and the strength of each component, we take as 
before the quantization axis z along the pair of Ce ions, and take the unit vector z. We obtain 
6;(M, M', z) oc 6 MM', since eq.(32) requires M' — M = m and Yijn{z) vanishes unless m = 0. 
This in turn gives the diagonal property: Bi{ya, v' a' ^ z) oc ?>vv'^ccj' ■ The projected interaction 
between a pair along the z axis is then given by 

D^^^-'\a,a' ,zR) = ^ S^z^a, i/a, i)i?;,(i/V, i/V, i)Q„,(fe,.i?). (55) 

1,1' 

For the nearest-neighbor interaction in the cubic lattice, the above choice of the quantization 
axis is consistent with the CEF states. On the other hand, the next-nearest-neighbor interac- 
tion in the Fg subspace requires the change of the quantization axis since the orbital index v 
is given for the cubic axis.^°) 

As was shown in ref,^°^ the high-symmetry around the z axis makes D^''y\a,a' , zR) 
independent of a, i.e., D^^'^'\a,a' ,zR) = D^^'^'\zR) . For any representation F, we can 
represent the relevant interaction matrix Dr{f^ , f^; zR) in terms of D^^'^'\ For example, 
components the of quadrupole-quadrupole interaction D^g is given by 

Dll{f,f;zR) =[D^+'+\zR) + D^-'-\zR)]l^, (56) 
Dl^gif, f; zR) =[D(+'-\zR) + D(-'+\zR)]/8. (57) 

Components D^g{f^, zR) with a = x,y,z turns out to be the same as Df^^f^, /°; zR) in 
consistency with the results in Fig.l. 

The sign of the anisotropy is clearly seen by retaining 1,1' up to two in the partial wave 
decomposition. This approximation is reasonable because the anisotropic range functions 



17/30 



J. Phys. Soc. Jpn. 



Full Paper 



Table IV. Tho (-xpausiou cocITick-uls of iuloractious Dr(f'\ f'\ :/?) l)y D'^''''''\zR). 



r 


compoucul D^^ 






4ul 




1/8 





4u2 


± 


1/32 


3/32 




II 


1/8 





3g 


± 





1/8 






1/8 





5g 







1/8 


2u 







1/8 


5u 




3/32 


1/32 









1/8 



with 1,1' ^ are smaller than |Qoo(a^)| at small \x\. We then obtain the compo- 
nents of the T3 quadrupole-quadrupole interaction as fohows: 



Drg{f,f;zR)c:, 



DTAf,f-zR)^ 



327r4 

327r4 



QooikpR) + ( ^ ) Q22{kFR) 



QooikpR) 



^] Q22{kFR) 



(58) 



(59) 



In general the hybridization of conduction electron causes fcrro-typc coupling between ions 
at short distance, and we have indeed Q22{x) < at x around 3. Then for the quadrupole- 
quadrupole interactions we have a relation Z?!^ < D^g at the nearest-neighbor distance. This 
is indeed satisfied by the results shown in Fig.l. 

Table IV shows interactions between operators within F by the coefficients of D^^''^'^ 
listed in the top column. Components D^°'{f^ , f^; zR) with a = x,y are shown with _L, and 
-Op(/°, f; zR) with II . From Table IV, it turns out that 1)^2 < ^ix2> and D^^ < D^. These 
two relations are consistent with the results in Figs 3 and 4. 



3.5 Anisotropy of the a -a Interaction 

It has been found that the exchange interaction between a pair of pseudo-spins a becomes 
isotropic in the case where the intermediate states of hybridization is either 4/° or 4/i4.i0'33) 
The isotropy follows independently of the local symmetry around the pair. As described in 
the main text, however, the next-nearest-neighbor exchange interaction in the Fg subspace is 
anisotropic. Here we clarify the condition under which the effective exchange becomes isotropic 
or anisotropic. 
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3.5.1 No orbital degeneracy 

Let the CEF ground state be F which has the Kramers degeneracy, but no orbital one. 
Within this CEF level, the on-site exchange interaction at a 4/^ site Ri can be written in the 
form 

i^exch (^) = ^ Xl^, (z)4^, {i)cva {i) , (60) 

aa' 

where X^^, (i) denotes either the transition from a' to a within the Kramers doublet, or the 
projection to a in the case of cr = a'. The conduction states are given in terms of the Wannier 
states. We note that this form of the exchange results only if the intermediate state is isotropic. 
The intersite interaction between Ri and Rj is determined by the average: 

{4aii)^ra'ij)), (61) 

and its Hermitian conjugate. With the time reversal symmetry preserved in the paramagnetic 
state the average is nonzero only if a = a' . Then the intersite interaction is proportional to 
the permutation operator 

Pij = ^X^^,(i)X^,^{j). (62) 

As in the main text wc introduce the pseudo-spin operator (Jx{i) = X^^ii) + X^^{i) and other 
components (Jy{i) and (7z(i) in an analogous way. Here the first component of the Kramers 
pair is indicated by j and the second one by |. Because of the identity 

^ii = ^k(0-o-(i) + l], (63) 

the effective exchange is found to be isotropic. We emphasize that the isotropic exchange is 
independent of the spatial symmetry around the pair, since it is the time-reversal symmetry 
that protects the isotropy. 

3.5.2 Orbital degeneracy 

Now we turn to the case of orbital degeneracy such as Fg. The exchange interaction now 
has orbital components specified by a, b. Provided the intermediate state of hybridization is 
isotropic, only the permutation of Fg components happens between the conduction states. 
Namely we obtain 

-ffexch(i) = JrYl ^aa 6a' (Oia' (O^raa (i), (64) 
aa'ab 

in the Wannier representation for conduction states. We now examine the possibility 

(4aT«cm(i))^0, (65) 

with a 7^ 6. In the case of the nearest-neighbor pair, the four-fold symmetry around the 
pair excludes the possibility, as discussed in §3.4. Then the F4U1 intersite interaction becomes 
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isotropic in the spin space. This situation was already noticed in ref.lO. 

However, if the second site measured from the first has the direction R other than the 
four-fold axis, we may have 

YlmiR) + 

with m = ±1 and ±3 for |m| < Z < 6 in eq.(32).^^) This is indeed the case of a next-nearest- 
neighbor pair where the first site is at the origin and the second one is at (a, 0, a). Prom eqs.(32) 
and (33), it is seen that a conduction electron including the component jz = ~l/2 from the 
first site can propagate to another state with = 5/2 in the second site with m = — 3, and 
to jz = —3/2 with m = 1. This situation realizes the possibility eq.(65). Note that eq.(65) 
does not necessarily mean the breakdown of the time reversal. The operators Tyax and Tyay, 
which are even under time reversal, describe the simultaneous flips of two kinds of quasi spins; 
orbital and Kramers pairs. 

Let us analyze a part of the exchange interaction which involves only the pseudo-spin 
component cr, but not the other component r. In the cubic symmetry this corresponds to the 
representation r4ui. The absence of the r component means that on-site process is diagonal 
with respect to the orbital index. However the Kramers index may have such combination 
with a 7^ a as 

a ± 

and its Hermitian conjugate. Expansion of the above form gives rise to an anisotropic term: 

a^i)-a-{j)-ayii).ayij), (67) 

in the pseudo-spin space cr. We understand in this way the results in Fig.2 which indeed show 
anisotropy in the pseudo-spin cr. 

We now discuss why there is no off'-diagonal component in Fig. 2(a). It is convenient to 
use the bond-density defined by eqs.(53) and (54) with the properties 

bi{M,M',R)^{-iy^-^''y2(^' ^ ^']Ylm{R), (68) 

m -M'J 

bi{-M',-M,R) = bUM,M',R), (69) 

Bi{ua, u'a', R) = Bl{v'a', va, R), (70) 
with a = —a. Then D^^^if^, /°, R) is proportional to 

Y,Qll'{kFR)Y.^ia')afs{an-ysBiiua,u'S,R)Bi,iu'j,u/3,R). (71) 

1,1' v,v' a/37<5 

Since R is in the x^;-plane, Yim{R^ is real with the azimuthal angle ^ = 0. Then using eq.(70) 
we see that Bi is also real, and hence eq.(71) is zero with Qn'ikpR) = QinikpR)- Therefore 
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^iuiif^^ f^i ^) vanishes when R is in the xz-plane. 
4. Summary 



We have derived the multipolar interaction between Ce ions from the Anderson-type model 
with orbital degeneracy using the free-electron model for the conduction electrons, with atten- 
tion to the electronic structures of the phases II and III in CeBg. As for intermediate states, 
we have considered not only 4/° but 4/^ Hund's-rule ground states, and have shown that the 
anisotropic exchanges which are expressible in terms of pseudo-spin operators arise. Prom the 
effective Hamiltonian, we have calculated coupling strengths for dipole-dipole, quadrupole- 
quadrupole, and octupole-octupole interactions. In this paper we have assumed the single 
spherical Fermi surface, which is evidently oversimplified if we compare with the experimen- 
tal situation in CeBg. Moreover there should be the Coulombic exchange in addition to the 
hybridization as another source of the intersite interaction. In spite of our simplification, the 
derived coupling constants give fairly consistent account of the AFQ order of Tr,g type in 
the phase II. The pseudo-dipole coupling constant of ^4^2 type and the anti-pseudo-dipole 
coupling constant of T^u type are both in consistency with the phenomenological model,^) 
provided 4/^ intermediate states are taken into account. 

In relation to the previous work,^°) we calculated the explicit spatial dependence of mul- 
tipolar couplings within free-electron approximation. The all coupling constants of multipolar 
interaction in this system have the same order of magnitude within this simplified model. Es- 
pecially in considering only 4/° intermediate state, our calculation is consistent with the pre- 
vious symmetry analysis on the couplings between nearest neighbors. The octupole-octupole 
interaction between type moments is as strong as the quadrupole-quadrupole interac- 
tion between T^g type moments. For next-nearest-neighbors, multipolar couplings become 
anisotropic even if only 4/° intermediate state is considered. 

As for the 4/^ intermediate states, we have considered only the Hund's-rule ground states. 
Inclusion of other states in the 4/^ multiplets should have a tendency to enhance the simul- 
taneous exchange of spin and orbit. This can be seen from the results in Table I. Namely, 
simultaneous spin and orbital exchange has additive contribution from both intermediate 
states with S = and S = 1. After projection to the Fg states, the orbital exchange is de- 
scribed by the pseudo-spin r, and the spin exchange by a. Although the spin-orbit interaction 
mixes the original orbital and spin quantum numbers considerably, the sign of the on-site ex- 
change is mainly determined by the antisymmetry of 4/^ states, and is robust against such 
mixing. After projection to the Fg states, the simultaneous spin-orbital exchange takes the 
form of Fsp, F4„2 and Fs^ on-site interactions. It is suggestive that these interactions seem 
relevant to realizing the phase diagram of actual CeBg.^^ 
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Appendix: Integrals involving spherical Bessel functions 

In this appendix we derive analytic expression of the range functions. We need to calculate 
integrals such as 

q(li,l2;z) = [ x"^ dx{ji^ix)(pi2{x) + ji2ix)(pi^{x)}, (A-l) 



JO 

where li and I2 are even integers, and ^^(x) is defined by 

'^'(^)^ / -2 j2Mx')dx'. (A.2) 

Jo - x/ 

The function Qi^i^ikpR) used in the main text is related to q{li,l2; kpR) by 

Qhi.ikpR) = {kFR)-\{hM\kFR). (A-3) 

Larsen has already calculated Qiii2{x) in the case of li = /2)'^^'* but we need Qi^i^ix) also with 
^1 / h. 

Using the result in ref.31, (piix) is given by 

m=l 



( \ ^ 



— (/ + m-l)!! 1 

xnnx) + > — -71 

^ ' 4-; (l-m)\\ X" 

m=2,4,... ^ ' 



(A-4) 



where ni{x) is the spherical Neumann function of order I. Spherical Bessel function of order I 
can be written by using trigonometric functions: 

ii{x) = Ci{x) cos X + Si{x)sm.x, (A-5) 

where Ci{x) and Si{x) are the polynomial expressions of x"^, while spherical Neumann func- 
tion of the same order is given by 

ni{x) = Ci{x) smx — Si{x) cos X. (A-6) 

We obtain the analytic expression for q{li,l2,z) by using eqs.(A-4), (A-5) and (A-6) and by 
performing the following partial integrals successively: 

sin X sin x If cos x 

ax = —7 — r H / ax j, (A-7) 

xP (jp - \)xP-^ p-l J xP-^ ^ ' 

cos X cos X 1 /" , sin x 
ax = —7 — T- / dx (A-8) 

The final results are tabulated in Table A-l. The functions in Table A-l are given in terms 
of the coefficients of the trigonometric functions and their related functions listed in the top 
row. 

Among the various range functions, Qoo{kFR) gives the classical RKKY interaction,^^) 

QooikpR) = 2TrF{2kFR), (A-9) 



/ 
/ 
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where 

F{x) = (xcosx — smx)/x'^. (A-10) 
Prom Table A-1, the asymptotic form of Qi^i^ikpR) for kpR — oo is obtained as:^) 

Qi.i.ikpR) {-lf^+'^y^Qoo{kFR). (A-ll) 
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Table A-1. Dimensionless range functions 2g(Zi, Z2; a;)/7r. The functions are given in terms of trigono- 
metric functions, Si(a;), and Si(2a;) (Si(a;) = dtsint/t). 



Table A-1. a The coefficient of a; "cos a:. 



h 


I2 X ^ cos(a;) 


X ' cos(a;) 


X cos(a;) 


X ^ cos(x) 


X ^ cos (a;) 









2 
4 
6 






-3465/4 


-105/2 
1575/8 


2 
2 
2 


2 
4 
6 




10395 


315/2 
2205/8 


9 

45/2 
-1071/8 


4 
4 
6 


4 
6 

6 21611205 


1091475/4 
-4584195/4 


3675 
-23625/8 
-255339/4 


175 
-11025/8 
-7119/2 


-5 
315/8 
189/4 



Table A-l.b The coefficient of x "sin a;, Si(a;) and Si(2a;). 



h 


^^2. X sin(a;) 


X * sin(a;) 


X ** sin(a;) 


X * sin(a;) 


X ^ sin(a;) 


Si(i;) 


Si(2i;) 









2 
4 
6 






-10395/4 


-105/2 
1575/8 


3 
-45 
210 


-3 
45 
-210 


2 
2 
2 


2 
4 
6 




-10395 


-315/2 
25515/8 


-9 

30 

-945/4 


3 
18 
-150 


-3 

-18 
150 


4 
4 
6 


4 
6 

6 -21611205 


-1091475/4 
33399135/4 


-3675 
751275/8 
324135/2 


1050 
12915/2 
5103/2 


115 
105/8 
-504 


10 
45 
21 


-10 
-45 
-21 



Table A-1. c The coefficient of a; "cos(2a;). 





li 


X ^ cos(2a;) a; ' cos(2a;) 


X ° cos(2a;) x' 


cos(2a;) x ^ 


cos(2a;) 


X cos(2a;) 
















1/2 





2 










-1/2 





4 










1/2 





6 






-3465/2 




-1/2 


2 


2 








-9 


1/2 


2 


4 






-315/2 


30 


-1/2 


2 


6 




-10395 


7245/2 


-63 


1/2 


4 


4 




-3675 


3325/2 


-100 


1/2 


4 


6 


-1091475/4 


278775/2 


-24885/2 


210 


-1/2 


6 


6 


-21611205 47806605/4 


-2819313/2 


95571/2 


-441 


1/2 


Table A-1. d The coefficient of a;"" sin(2a;). 


h 


I2 


a:-^"sin(2a;) a:"** sm(2a:) 


X sin(2a;) 


X *sin(2a;) x 


sm(2a:) 


sin(2a::) 
















-1/4 





2 










7/4 





4 








105/2 


-21/4 





6 






10395/4 


-630 


43/4 


2 


2 








9/2 


-13/4 


2 


4 






315/4 


-120 


27/4 


2 


6 




10395/2 


-34965/4 


945 


-49/4 


4 


4 




3675/2 


-13125/4 


505 


-41/4 


4 


6 


1091475/8 


-1006425/4 


202545/4 


-2100 


63/4 


6 


6 


21611205/2 -163066365/8 


19322415/4 


-1210167/4 


10899/2 


-85/4 
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(a) 0.002 




Fig. 1. Spatial dependences of quadrupole-quadrupole interactions D^g and D^g with R = (0, 0, R). 
Results with different intermediate states are shown as follows: (a) intermediate states on both 
sites; (b) intermediate states on one site and on the other site; (c) intermediate states 
on both sites. Arrows show kpRnn- 
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Fig. 2. The spatial dependence of dipole-dipole interaction -D4„i. The part with kpR < 4 shows the 
case of R = {0,0, R), and the other part with kpR > 4 gives the result for R = {R,0,R)/^/2. 
Arrows show kpRnn and kpRnnn- 
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Fig. 3. The spatial dependence of dipole-dipole interaction -D4„2- The cases (a), (b), (c) are the same 
as in Fig.l. And the defintion of R are the same as in Fig. 2. Arrows show kpRnn and kpRnnn- 
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Fig. 4. The spatial dependence of octupole-octupole interactions D2u and D^^. The cases (a), (b), 
(c) are the same as in Fig.l. And the definition of R are the same as in Fig.2. Arrows show kpRnn 

and kpRnnn- 
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